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Gauge-Invariant Initial Conditions and Early Time Perturbations in Quintessence 
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We present a systematic treatment of the initial conditions and evolution of cosmological per- 
turbations in a universe containing photons, baryons, neutrinos, cold dark matter, and a scalar 
quintessence field. By formulating the evolution in terms of a differential equation involving a ma- 
trix acting on a vector comprised of the perturbation variables, we can use the familiar language 
of eigenvalues and eigenvectors. As the largest eigenvalue of the evolution matrix is fourfold de- 
generate, it follows that there are four dominant modes with non-diverging gravitational potential 
at early times, corresponding to adiabatic, cold dark matter isocurvature, baryon isocurvature and 
neutrino isocurvature perturbations. We conclude that quintessence does not lead to an additional 
independent mode. 

PACS numbers: 98.80.Bp, 98.80.Cq 



I. INTRODUCTION 

The advent of high precision data of the cosmic 
microwave background (CMB) anisotropics permits de- 
tailed tests of the composition and shape of the pri- 
mordial density fluctuations. The most popular models 
of inflationary cosmology predict adiabatic fluctuations 
HI El El More elaborate models lead to an admixture of 
adiabatic and isocurvature fluctuations [a, |6| . The time 
evolution of adiabatic and non-adiabatic fluctuations is 
well understood for a universe composed of radiation, 
baryons, cold dark matter (CDM) and neutrinos [7|. In 
the context of quintessence 0, 0, , the behaviour of 
the field fluctuation has been studied in several works 
fill Hp. Ilii . FLU Il5l | . Initial conditions have been proposed 
in |16| for the case of negligible quintessence contribution 
in the early universe. We present here a systematic treat- 
ment of initial conditions for quintessence models which 
differs from that of |l6j | in approach and interpretation. 

Our basic setting assumes that small deviations from 
homogeneity are generated during a very early stage of 
the big bang, typically an inflationary epoch. During 
the following radiation dominated period the wavelength 
of the relevant fluctuations is far outside the horizon. 
Apart from this, we will not use any further constraint 
on the primordial fluctuations. Only the spectra of a cer- 
tain number of "dominant" modes can possibly influence 
events such as emission of the CMB and its anisotropies 
since the other modes decay. The information about 
these dominant modes therefore constitutes the initial 
conditions for practical purposes. Primordial informa- 
tion beyond the dominant modes is effectively lost and 
not observable. The detailed time of specification of the 
initial conditions is therefore irrelevant as long as it is 
much shorter than the time of matter-radiation equality. 

During the period relevant for the discussion of the 
initial conditions the universe is radiation dominated. 
However, our approach allows for the presence of scalar 



fields which evolve like radiation at early times or are 
subdominant. Consequently, our results hold for a wide 
class of quintessence models, including those with non- 
negligible Cl q at early times 0. In fact, we only use 
a "tracking" property for the background of ho- 
mogenous quintessence, namely that its equation of state 
w q = Pq/ p q is almost constant and determined only by 
the energy densities of the radiation and matter compo- 
nents. The parameters w q and Q q = I — Q m — Q v — f2 7 
will therefore be the only parameters of the quintessence 
model that influence the early time evolution of small 
fluctuations. This makes our analysis model independent 
to a large extent. 

We will formulate the evolution equations for the per- 
turbation variables as a first order differential matrix 
equation: 



d In. 



U = A(x)U, 



(1) 



where the vector U contains all perturbation variables 
and the matrix A(x) encodes the evolution equations. In 
doing so, we relate the problem of finding initial condi- 
tions and dominant modes to the familiar language of 
eigenvalues and eigenvectors. This formulation makes 
"mode-accounting" transparent by counting the degen- 
eracy of the largest eigenvalue. We find four dominant 
modes that remain regular at early times. For phys- 
ical reasons, we choose a basis using adiabatic, CDM 
isocurvature, baryon isocurvature and neutrino isocur- 
vature initial conditions. As we will show, adiabaticity 
between CDM, baryons and photons implies adiabaticity 
of quintessence. There is therefore no pure quintessence 
isocurvature mode. In addition, using the matrix formu- 
lation reveals facets of the modes that otherwise remain 
obscured. 

In order to avoid the appearance of gauge modes, we 
will use the gauge-invariant formalism [19l 1201 l2ll |22|. 
In contrast to earlier work, we find it more appropri- 
ate to specify the initial conditions and time evolution 
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of the quintessence field in terms of the gauge-invariant 
density contrast and velocity, thus unifying the language 
for all species. As anticipated, the quintessence density 
perturbation remains constant at super-horizon scales for 
adiabatic initial conditions. In contrast to this, the field 
fluctuation follows a simple power law in conformal time 
that only depends on the quintessence equation of state. 

We will proceed as follows: in section [D] we give the 
gauge-invariant perturbation equations for a radiation- 
dominated universe containing radiation, cold dark mat- 
ter, neutrinos, baryons in the tight coupling limit and 
tracking quintessence. We express the evolution in ma- 
trix form in III Bl In section UlI Al we classify the modes 
and determine them in sections IlII Bl HITOl and lni Dl To 
illustrate the effect of non-adiabatic contributions to the 
CMB spectrum, we plot a few spectra for different initial 
conditions in section llVl A summary of our findings is 
given in section^ In Appendix^ we derive the pertur- 
bation equations used in detail, while Appendices iBl and 
10 discuss supplementary issues. 



II. THE PERTURBATION EQUATIONS 



In the following we adopt the gauge- invariant approach 
as devised by Bardeen [lS|. It is not difficult to obtain 
the initial conditions in any gauge from the corresponding 
gauge-invariant quantities given here. In Appendix A, we 
summarize the definitions of the perturbation variables 
and sketch the derivation of the evolution equations. It 
turns out that the evolution is best described as a func- 
tion of x = kr, where r is the conformal time and k the 
comoving wavenumber of the mode. We assume that at 
early times, the universe expands as if radiation dom- 
inated. This assumption is well justified for small fl q 
at early times, as well as for potentials that are essen- 
tially exponentials at the time of interest, regardless of 
£l q . The assumption is certainly not justified for mod- 
els in which quintessence is dominating the universe at 
early times with equation of state w q ^ 1/3. For such 
(slightly exotic) models, the following steps would need 
to be modified. 
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n:, 
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■=Vv ~ 2IL,, 
5 



A' = 



vl = 3fi„n„ 



3{Wg - 1) 

3 



A q + 3(1 + W q 



3K - 1) 
A„ 



(1 + W q )Vq 



Vq+W, 



(8) 
(9) 

(10) 
(11) 



1 + Wq 

with the gauge-invariant Newtonian potential "J given by 



E n a (A a + 3(i + Wa )v a ) 
E 3(i + w a )n a + if 

a=c,b,7 ,v.q 



tl„ILv. (12) 



We denote the derivative d/dlnx with a prime. The 
gauge-invariant energy density contrasts A a , the veloc- 
ities V a and the shear 11^ are the ones found in the lit- 
erature 0, |2jJ, H3, except that we factor out powers 
of x from the velocity and shear defining V = V/x and 
Ii v = x~ 2 Ilv. This factoring out leads to the particu- 
larly simple form of the system of equations for x <C 1 
(see also Appendix lA"|l . It does, however, exclude modes 
with diverging \& at early times such as a neutrino veloc- 
ity mode |23|. The index a runs over the five species in 
our equations, namely cold dark matter, baryons, pho- 
tons, neutrinos and quintessence, denoted with the sub- 
script q. We assume tight coupling between photons and 
baryons. The equation of state w = p/p takes on the 
values w c = Wb = 0, n? 7 = w v = 1/3 and w q is left as 
a free parameter. Equations J5J, © and J7J) can 
be regarded as continuity relations between the density 
fluctuations and the velocity. We obtain equations 11U|) 
and HI IB from the perturbed Klein-Gordon equation of 
the quintessence scalar field expressed in terms of A q 
and V q , the energy density and velocity perturbations as 
defined in Appendix lAl 



A. Full Set of Equations 

Assuming tracking quintessence we obtain the follow- 
ing set of equations (for a derivation, see Appendix |A"|> : 



A^ - -x 2 K, (2) 

V' c = -2V C + V, (3) 

a; = ~x%, (4) 

% - jA 7 - % + n v ii u + 2*, (5) 

A^ = -x 2 %, (6) 

K = ~^Vu, (7) 



B. Matrix Formulation and Dominant Modes 

Conceptually, it is convenient to note that the above 
set of equations can be concisely written in matrix form 
according to Equation (JIJ where the perturbation vector 
is defined as 

U T = (A c , V c , A 7 , %, A b , A,, V v , IL V , A q , V q ). (13) 

The matrix A(x) can easily be read off from equations 
(|2*|)- (|llfl . This enables us to discuss the problem of spec- 
ifying initial conditions in a systematic way. 

The initial conditions are specified for modes well out- 
side the horizon, i.e. 1 « 1. In this case, the r.h.s. of 
equations J2J > ® > ® and (7J can be neglected, provided 
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V a does not diverge oc x~ 2 or faster for x 2 — > 0. The 
evolution matrix A(x) loses any explicit x dependence 
for x 2 — > 0. Yet, it still depends on x via terms involving 
fl c , fib and fl q . By our assumptions on quintessence, the 
term involving fl q is either a constant (for w q = 1/3) or 
negligible (yet, in Appendix El we extend the treatment 
to include models with considerable fl q and w q =/= 1/3.) 
In both cases fl q can be approximated by a constant 
(fl q = for w q < 1/3) and fl c , fib vanish oc x. In leading 
order, the matrix A becomes therefore x-independent for 
very early times. In fact, the general solution to Equa- 
tion in the (ideal) case of a truly constant A would 
be 

^) = E c *(r) A ' E7(i) ' (14) 

i 

where U^> are the eigenvectors of A with eigenvalue A; 
and the time independent coefficients c; specify the ini- 
tial contribution of towards a general perturbation 
U. As time progresses, components corresponding to the 
largest eigenvalues A^ will dominate. Compared to these 
"dominant" modes, initial contributions in the direction 
of eigenvectors U^ 1 ' with smaller Re(A,) decay. It there- 
fore suffices to specify the initial contribution c; for the 
dominant modes, if one is not interested in very early 
time behaviour shortly after inflation. In our case, the 
characteristic polynomial of A(x) indeed has a fourfold 
degenerate eigenvalue A = in the limit x 2 — > 0, inde- 
pendent of fl c , fib and flq. 1 While it is not feasible to 
obtain the remaining six eigenvalues by analytic means, 
we have checked numerically for a wide range of f2 7 , fl Vl 
fib, fl c , fl q and w q that the remaining eigenvalues have 
indeed negative real parts and contributions from the 
corresponding eigenvectors towards a general perturba- 
tion U will therefore decay according to Equation (|14f> . 
We can improve the analytic description of the dominant 
modes by taking corrections oc x into account. 

As fl c oc fib oc x, it is appropriate to split A(x) accord- 
ing to the scaling with x, 

A = A + xAi, (15) 

where Aq and A\ are constant and x A\ contains the 
small, time-dependent corrections from terms involving 
fl c and fib- We may also write 2 the eigenvectors as a 
series in x, 

U = U + xU 1 . (16) 
Inserting Equations l(T5 )l -l|TS )l in Equation JIJ, we get 

A) U = 0, (17) 



1 For w q = 1 we find another eigenvalue with A = 0. We will 
ignore this special case in what follows. 

2 This form is not an ansatz, but dictated by Equation ft) , once 
the dependence of A(x) on x is given. 



and 

U x = -(A - l^AxUo. (18) 

Equation (|18|l is easy to solve, once Uq has been deter- 
mined (we discuss the possibility of a vanishing Uq in 
Appendix [UJ . We see from eqn. H17|) that to constant 
order the solutions of eqn. are indeed given by eigen- 
vectors to the eigenvalue A = 0. We should emphasize 
that the vectors Uo do not evolve in time if their corre- 
sponding eigenvalues are A = 0. Thus, the perturbations 
remain constant in the super-horizon regime during ra- 
diation domination in this approximation. If we include 
the next-to-leading order contribution to U, the eigen- 
vectors do evolve and we can no longer apply eq. l|14fl . 
These corrections are, however, small as long as we are 
deep in the radiation dominated era due to the small con- 
tributions of baryons, radiation and quintessence during 
this era. Given a set of initial conditions in the form of 
coefficients for the four dominating modes at Zinitiai we 
can find the perturbations at some later time (provided 
the modes are still super-horizon sized and we have radi- 
ation domination). In leading order, the coefficients will 
remain the same while in next-to-leading order we can 
use the evolution of U to compute the coefficients for 
z < Zi n itiai- If initial conditions are specified with accu- 
racy of next-to-leading order one therefore has to specify 
Zinitiai as well. In leading order this is unecessary for z 
in a wide range long before last scattering. 

C. Constraint Equations to Leading Order 

Equation l|17|) is equivalent to setting the l.h.s. of 
Equations (|2*|)- (|ll|l equal to zero and using fl c — fib = 
x 2 = 0. Then Equations 0, igj, © and (0 are auto- 
matically satisfied (provided V a does not diverge oc x~ 2 
or faster), and Equations ll3l.(|51).llS|l- Hll|) yield non-trivial 
constraints for the components of Uq : 





2V C 


- * 


= o, 


(19) 


1/4A 7 - 


% + fi v ii v - 




= o, 


(20) 


1/4A„ - 


V V + flyliy ~ 


-2* 


= o, 


(21) 




8/5K - 


2n, 


= o, 


(22) 


2>fl v Tl v + A q /(1 - 


h»,) +3Vq - 


-3* 


= o, 


(23) 


3fi„n„ + A,/(l 


+ Wq) + Vq - 


-4* 


= 0. 


(24) 



In the above, all quantities are considered only to con- 
stant order, (we have omitted the subscript '0' for no- 
tational convenience.) In particular, there is no contri- 
bution of CDM and baryons to ^ at constant order. 
Note that, apart from w q , no model-specific parame- 
ters occur in any of these equations so the modes will 
be independent of the type of quintessence as long as 
the scalar field is in a regime with approximately con- 
stant Wq. We note that for w q substantially smaller than 
1/3 the quintessence fraction fl q changes with time. By 
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the assumption that the universe expands as if radiation 
dominated, the quintessence contribution would however 
be small in this case and its contribution to '5 can be 
neglected (see Appendix IU1 for an extended discussion). 

We mention that for w q = 1/3, quintessence evolves 
the same way as radiation, therefore Q q does not change 
in this case. If w q = —1/3, quintessence has the same 
influence on the scale factor a as a curvature term in an 
open universe. However, the geometry is still flat and one 
can distinguish an open universe from this quintessence 
model by measuring the position of the first acoustic peak 
in the CMB. 



III. THE MODES IN DETAIL 

A. Classifying the modes 

While any basis for the subspace spanned by the eigen- 
vectors with eigenvalue A = can be used to specify 
the initial conditions, it is still worthwhile to use a basis 
that is physically meaningful. Following the existing lit- 
erature, we use the gauge-invariant entropy perturbation 

M 



s, 



A 



a:/3 



1 + W a 1 + Wp ' 



(25) 



between two species a and /?, as well as the gauge- 
invariant curvature perturbation on hyper-surfaces of 
uniform energy density of species a 0, la.123. |2j| 



5p a 



3(1 + w a )p a 



(26) 



in order to classify the physical modes. On slices of uni- 
form total energy density, the curvature perturbation is 
correspondingly 



Ctot = H 



1 



Ea 3 ( 1 + W "W 



(27) 



In our variables, these expressions take on the manifestly 
gauge-invariant form 



Co 



A, 



3(1 + IBa) ' 



J2 a 3(1 + W a )£l a 



(28) 



If Ctot = 0, energy density perturbations do not generate 
curvature. It is therefore clear that such a perturbation is 
a perturbation in the local equation of state. One should 
note that the definition of Ctot is different from that of 
El: 



B. The Adiabatic Mode 

The first (rather intuitive) perturbations one would try 
to find are adiabatic perturbations, which are specified by 
the adiabaticity conditions S a: p = for all pairs of com- 
ponents. In our case, this results in eleven constraints 3 
for the ten components of Uq. It is a priori not clear that 
this has a solution so we will not include quintessence in 
the adiabaticity requirement. Requiring adiabaticity be- 
tween CDM, baryons, neutrinos and radiation, 



A„ = A 7 = -A c = -A h , 



(30) 



and using the six constraint Equations ijT§|) - (|?lJl . we ob- 
tain 



/ A c \ 






( 


3/4 




Vc 






(- 


■5/4) V 




A 7 








1 




V-f 






(- 


■5/4) V 




A b 
A, 




= C 




3/4 
1 




Vu 






(- 


■5/4) V 




flu 








-V 










3(1 


+ w q )/A 




V v q J 


adiabatic 




V (- 


-5/4) V 


) 



(31) 



where V = (15 + 451^)^ 1 and C is an arbitrary con- 
stant. From A g /A 7 = 3(1 + w g )/4 we conclude that 
quintessence is automatically adiabatic if CDM, baryons, 
neutrinos and radiation are adiabatic, independent of the 
quintessence model for as long as we are in the tracking 
regime. As all components are non-vanishing, we do not 
quote the next to leading order contributions from x U\ . 



C. Neutrino Isocurvature 

Having found the adiabatic vector, one could specify 
three additional linearly independent vectors satisfying 
the constraint Equations (|19l) - (|24|l . This would com- 
plete the basis. It is, however, appropriate to choose 
modes that may be generated by physical processes. 
These modes are in general not orthogonal but span the 
eigenspace of A — 0. Modes that may be generated by 
physical processes are isocurvature modes. A given mode 
is an isocurvature mode, if the gauge-invariant curvature 
perturbation Ctot vanishes, i.e. Ctot = 0. In order to dis- 
tinguish different isocurvature modes from one another, 
we require that the other species are adiabatic with re- 
spect to each other, i. e. S a: p = except for quintessence 
and one species er, which has non- vanishing Sa-^. 



Cmfb 



2H- 1 i> 
3" 



(1 



(29) 



However, one may verify that this quantity coincides with 
Ctot in the super-horizon limit for a flat universe |26| . 



3 Without requiring quintessence to be adiabatic, we have six con- 
straints from equations 1191 - 1241 plus three constraints from eq. 
OHJ Plus one constraint from the overall normalization, which is 
fixed by choosing a specific value for A 7 . 
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Let us first consider the neutrino isocurvature mode. 
For this, we require that CDM, baryons and radiation are 
adiabatic, while S u: ~ 7^ and that the gauge-invariant 
curvature perturbation vanishes: 



C tot = 0, A c = A b = -A 7 . 
Using this and Equations (|19[l - i|24[l leads to 



(32) 



/ A c \ 






( 


3/4 \ 


Vc 










A 7 








1 


% 








(a, + n„ + f ) v 


A, 




= C 




3/4 


A, 








v u 


















-svsi 7 /n v 













V % ) 


neutrino iso. 




V 





It is important to note that we did not 



(33) 
require 



quintessence to be adiabatic. One can see from the neu- 
trino isocurvature vector that A q = 0, and as a con- 
sequence quintessence is not adiabatic with respect to 
either neutrinos, radiation, baryons or CDM. Hence, we 
could just as well have labeled this vector "quintessence 
isocurvature" . We cannot require adiabaticity between 
neutrinos, CDM, baryons and radiation and hope to ob- 
tain a "pure" quintessence isocurvature vector since, as 
we have seen in the discussion of the adiabatic mode, 
these requirements lead to quintessence being adiabatic 
as well. 



D. CDM isocurvature and baryon isocurvature 

The CDM isocurvature mode is characterized by 
Sc-.-y 7^ 0, C,tot — and adiabaticity between photons, 
neutrinos and baryons: 



Q ot = 0, A 7 = A„ = -A b . 



(34) 



Using this and Equations yields 

£/ T (CDM iso.) = (1,0,0,0,0,0,0,0,0,0). (35) 

This vector fulfills (tot = + 0(£l c ), which is in line 
with our approximation since fi c <C 1. Similarly, for the 
baryon isocurvature mode, we require S b:l 7^ 0, Ctot = 
and adiabicity between photons, neutrinos and baryons. 
The resulting vector reads 

[/^(baryon iso.) = (0,0,0,0,1,0,0,0,0,0). (36) 

As all but one of the components of Uo are vanish- 
ing for CDM isocurvature and baryon isocurvature, we 



use Equation 1|18|) to obtain the next to constant order 
solution for CDM isocurvature 



/ A c \ 

Vc 

A 7 

V, 
A 6 
A„ 
V v 
ft 

V v q J 



c 



CDM iso. 



/ 1 \ 

n c (4fi„ - 15)W/12 



-(15/4)fi c W 



-(15/4)fi c W 
-2tt c U 
O c (15 + 20^(1 + w q )U 

V n c uv 



where U = (30 + 4O,,)- 1 

40,(3w g -l)]/[36K-l)] 

isocurvature 



/ A c \ 

V c 
A 7 
V 



and V 



[105 - 45w q 



J 

(37) 



"7 
A 6 
A, 

Vu 

A„ 



C 



aryon iso. 



Similarly, we find for baryon 



\ 

n 6 (4f2 I/ -15)W/12 



-(15/4)n b W 

1 



-(15/4)fi b W 

-2n b u 

n b (15 + 2Slu)(l + w q )U 

n b uv 



J 

(38) 

Note that these vectors are not constant since £l b and 
il c both evolve in time. We observe that the corrections 
to U are indeed proportional to fl c or fl b as expected. 
This result holds for all tracking quintessence models 
with w q = 1/3 or w q < during the radiation domi- 
nated period. For intermediate values < w q < 1/3 
the deviation from the leading behaviour scales cx x a , 
a < 1. Obviously, the adiabatic, CDM isocurvature, 
baryon isocurvature and neutrino isocurvature vectors 
Uo are linearly independent. We have therefore iden- 
tified four modes corresponding to the fourfold degener- 
ate eigenvalue zero of A(x). These four vectors span the 
subspace of dominant modes in the super- horizon limit, 
and there are no more linearly independent vectors that 
satisfy the constraints l|19|) - 124fl . Arbitrary initial per- 
turbations may therefore be represented by projecting a 
perturbation vector U at initial time into the subspace 
spanned by the four aforementioned vectors, as this is 
the part of the initial perturbations which will dominate 
as time progresses. 

Figure ^ demonstrates that the early time behaviour 
is well described by our analytic formulae. The analytic 
results agree very well with the simulation for early times, 
when the mode is outside the horizon. In the lower graph, 
we plot the equation of state w q . The quintessence model 
used is parameterized by an equation of state w q (a) = 
—0.95 + 0.75(1 — a), leading to w q (early) = —0.2 and 
according to (|B3|) . X cx r 8 . This differs from reference 
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FIG. 1: Gauge-invariant energy density perturbation A q and 
quintessence field fluctuation X as simulated (straight and 
dashed-dotted lines), compared to the analytic solution of 
Equations I3H and <B3t (dashed and dotted lines) as a func- 
tion of conformal time r for adiabatic initial conditions. Radi- 
ation and matter equality corresponds to r = 109Mpc. Shown 
is the mode for k — 0.1 Mpc -1 and the cosmological parame- 
ters have been fig/i 2 = 0.022, h = 0.7, = 0.3, = 0.7. 



M 4 

We see that including quintessence does not add a new 
dominant mode. The two additional modes added by 
the fluctuations of the scalar field are both subdomi- 
nant and decay with negative eigenvalue A;. This is 
due to the fact that none of the perturbation equa- 
tions for quintessence equate to zero in the superhorizon 
limit. This holds for non-tracking quintessence models as 
well. Let us investigate this in detail. For all the other 
fluid components, A^ = in the super-horizon limit, 



but for quintessence we get from eq. (| A29|) that A' q 
^ c s(?) — w q)^-q ~ 3w q T q . For tracking quintessence, we 



obtain from equation l|A46|) that cj/ g \ = w q and we find 



A; = -3w, 



(39) 



Since T q does not vanish except for w q — 1 (see eq. 
(|A45(l ). this does not equate to zero. 5 Hence, due to 
the non-vanishing entropy perturbation of quintessence 
there is no additional dominant mode. 6 



4 In it is stated that the quintessence fluctuation in Newtonian 
gauge scales oc r 2 for adiabatic initial conditions. This does not 
ag ree with our results in Appendix [E] Actually, equation (101) 
of includes a factor ip t o, which, interpreted as a dynamical 
quantity dip/dt (and not fixed at some initial time to), leads 
to a power law in r which is then consistent with our result of 
Appendix ITU 

5 Note that w q = does not lead to A' g = 0. 

6 We have not yet investigated the relationship between decaying 



IV. ISOCURVATURE INITIAL CONDITIONS 
AND THE CMB 

We illustrate the influence of different initial condi- 
tions on the CMB with an example. For an analysis 
of experimental data and a possible isocurvature contri- 
bution to the CMB we refer the reader to [27l I28L l29| . 
Here, we merely whish to show the qualitative features 
of the different modes. We use a modified version of CM- 
BEASY [30L l3l| to compute CMB spectra corresponding 
to different initial conditions for an early quintessence 
cosmology with parameters as in model A of |l7|. We 
set the spectral index of the isocurvature modes iden- 
tical to the spectral index of the pure adiabatic mode, 
n s = 0.99. The resulting spectra are plotted in Fig. 
The spectrum of the pure CDM isocurvature mode de- 
cays quickly when going to small scales as has been found 
in previous works [321 l33l l34j |. The neutrino isocurva- 
ture mode shows prominent peaks at higher multipolcs 
than the adiabatic mode with different peak ratios. For 
the mixed initial conditions with only small isocurvature 
contribution, the shape of the curve remains more or less 
the same. A small admixture of isocurvature fluctua- 
tions leads to a decrease of power at larger multipolcs if 
the overall normalization is fixed at I — 10. Comparison 
with the WMAP data in the same figure shows that non- 
adiabatic initial perturbations are strongly constrained. 
Clearly, pure isocurvature initial conditions arc inconsis- 
tent with CMB observations. 



V. CONCLUSION 

We have investigated perturbations in a radiation- 
dominated universe containing quintessence, CDM, neu- 
trinos, radiation and baryons in the tight coupling limit. 
The perturbation evolution has been expressed as a dif- 
ferential equation involving a matrix acting on a vector 
comprised of the perturbation variables. This formula- 
tion leads to a systematic determination of the initial con- 
ditions. In particular, we find that due to the presence 
of tracking scalar quintessence no additional dominant 
mode is introduced. This fact is beautifully transparent 
in the matrix language. Indeed, contributions of higher 
order in x = fcr towards a perturbation vector U can 
easily be determined by solving a simple matrix equa- 
tion once the constant part of U has been determined. 

In total, we find four dominant modes and choose 
them as adiabatic, CDM isocurvature, baryon isocurva- 
ture and neutrino isocurvature. For the neutrino isocur- 
vature mode, quintessence automatically is forced to 
non-adiabaticity. Hence, we could have as well labeled 
the neutrino isocurvature mode as quintessence isocur- 
vature. To demonstrate the influence on the cosmic mi- 



quintcssence modes and the background evolution. 
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FIG. 2: CMB Temperature spectra as a function of multipole / in an early quintessence cosmology. The pure adiabatic (AD), 
CDM isocurvature (CI), neutrino isocurvature (NI) mode and three different combinations of these dominant modes are plotted. 
For comparison with experimental data we also give the WMAP measurements of the CMB 0. The spectrum of the pure 
baryon isocurvature mode is essentially identical to that of the pure CDM isocurvature mode. All spectra have been normalized 
to the same power at / = 10. 



crowave background anisotropy spectrum, we have calcu- 
lated spectra for all modes. Clearly, non- adiabatic contri- 
butions are severely constrained by the data. A detailed 
study may provide ways to put additional constraints on 
quintessence models or tell us more about the initial per- 
turbations after inflation. 
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APPENDIX A: GAUGE-INVARIANT 
PERTURBATION EQUATIONS 

In this Appendix we will explain the derivation of equa- 
tions [t2")l- ljT2T) in detail. 



1. The general story 

First, we briefl y s ummarize the gauge- invariant ap- 
proach of |19t l20l [23] ■ Perturbing a homogenous Fried- 



man universe, one classifies fluctuations according to 
their transformation properties with respect to the ro- 
tation group. In flat spacetime, we may expand the per- 
turbation variables in terms of harmonic functions |35| . 
With Q t i = dQ/dx 1 one defines 

Q i {k,x) = -k- 1 Q{k,x), i (Al) 

and 

Q i0 {k,x) = k- 2 Q{k.x)^ + i<%Q(fc,a;), (A2) 

where the Q(k,x) are eigenfunctions of the Laplace- 
operator, V 2 Qk(x) = —k 2 Qk(x) and in spatially flat uni- 
verses Q = exp(ikx). As modes with different k decouple 
in linear theory, we will not display the fc-dependence of 
Q in the following. The scalar parts of vector and tensor 
fields can then be written as 

Bi = BQi (A3) 

and 

Hij =H L Q5 ij +H T Q ij . (A4) 

respectively. 

In this work, we are only interested in scalar fluctua- 
tions because scalar quintessence will not influence vec- 
tor or tensor modes. The most general line element for a 
perturbed Robertson- Walker metric may be written as 

ds 2 = a(r) 2 [-(l + 2A)dr 2 - 2B l drdx l + 

(Sij + 2H i3 )dx i dx j } 1 (A5) 
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where in the scalar case B t and Hij are given by equations 
(|A3I) and IIA4I) . The gauge transformation of a tensor T 
is given by [H HI Hi |s3 



f{x) = T(x) - L t f, 



(A6) 



where L e is the Lie derivative. The transformation vector 
e can be decomposed as 



T(t)Q(x), 



i l = x l + L{t)Q\x), 



(A7) 



(A8) 



where L and T are arbitrary functions of r. The transfor- 
mation properties of the metric perturbations are given 
by pHU 



4 = A-HT-T, 

B = B + L + kT, 

Hl = H L -HT- ~L, 

Ht = Ht + kL, 



(A9) 
(A10) 

(All) 

(A12) 



where a dot denotes the derivative with respect to confor- 
mal time r and H. = a(r)/a(T). The functions L and T 
can be used to eliminate two of the metric perturbations. 
Popular choices are A = B = for the synchronous gauge 
and B = Ht = for the longitudinal gauge. 

From equations ljA9(l - i|A12(l one can construct the 
gauge-invariant Bardeen potentials |19| 



* = A - Uk~ x a - kr x b, 
$> = H L + -H T - HArV, 



(A13) 
(A14) 



with a = k~ l H T - B. It is worthwhile to note that in 
longitudinal gauge, for which B — Ht = cr = 0, the 
perturbed metric takes on the simple form 



dsl long . ) =a(r) 2 [-(l + 2*Q)dT 2 



(l + 2<$>Q)5 ij dx i dx j ]. (A15) 



p = (8ttG) 1 / 2 denoting the reduced Planck 



With M P 
mass, Einstein's equation reads 



(A16) 



traceless part by 7T*j- = p LT Q z j. Therefore the compo- 
nents of the energy momentum tensor are 



r 



T 



= -p(l + 5Q), 

= -p(l + w)vQ\ 

T° = p(l+w)(v-B)Q i , 

T i j = p [(l + n L Q)S) + UQ)] 



(A18) 
(A19) 
(A20) 
(A21) 



Given the gauge- transformation properties of S, v and 
ttl |l9ll20ll2lll23.l30| . one can construct gauge- invariant 
quantities for the energy density contrast A, the velocity 
V and the entropy perturbation T. These are given by 



A 
V 

r 



S + 3(1 + w) [H L + ^H T 

v - k- l H T , 
c 2 

n L d. 

w 



(A22) 
(A23) 
(A24) 



Here, c 2 = dp/ dp is the adiabatic sound speed. From 
the conservation of the zero component of the energy 
momentum tensor V^Tq = we obtain 



Pn 

Pa 



-3(1 + w a )H, 



(A25) 



where w — p/p is the equation of state of the particu- 
lar species. The perturbed Einstein equations in gauge- 
invariant variables are [HL I20L I2I I2I l30| : 

a 2 pA = 2Mpk 2 $ 

-3a 2 p (1 + w) (Hk~ l V - $) , (A26) 

a 2 (p+p)V = 2M^fc(w-l), (A27) 



a 2 pn = -M 2 5 fc 2 ($ + *), 



(A28) 



In the above, it is understood that the quantities A, V 
and LT are the sum of the contributions of all species a. 
Using w = (c 2 — w)p/p and 1|A25|I we get from = 
that 

A + 3(c 2 - w)HA + kV(l + w) + 3HwT = 0, (A29) 



and from 







V = H(3c 2 s 
c 2 k 



1 



1)V + Jfc[* -3c 2 $ 

wk r 2 

a + r- -n 

l + w 3 



(A30) 



where the energy momentum tensor of a perfect fluid is 
given by 

T% = p&l + {p + p)u»u v + (A17) 

The covariant 4- velocity is itj = cl(v(t) — B) Qj. We de- 
fine the energy density contrast (5 by p = ,5(1 + 6(t))Q, 
the spatial trace by pS^ = p(r)(l + ttl{t)Q) 8^ and the 



2. Gauge-invariant quintessence perturbations 

The scalar quintessence field is decomposed into a 
background and fluctuation part according to ip(r, x) = 
<P(t) + x{t,x). The fluctuation can be promoted to a 
gauge-invariant quantity by defining the gauge-invariant 
quintessence field fluctuation X = x ~~ t pk~ 1 o. The field 
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Symbol 


Meaning 


Equation 


o 

^ "SpCClCS 


fraction of total energy density 


n.a. 


^ "species 


fraction of total energy density today 


n.a. 


a 


scale factor of the universe 


n.a. 


T 


conformal time: dr = At /a 


n.a. 


A' 


wavenumber ot mode 


n.a. 


X 


k'T 


n.a. 




derivative w.r.t conformal time 


n.a. 


l 


derivative w.r.t. x-\ — 

dx 


n.a. 


rt 


cl/(i 


n.a. 


A 


gauge- inv. density contrast (A g of [20]) 


llA"22Tl 


V 


gauge-invariant velocity 




n 


shear 


ilMl 


V 


reduced velocity: V = x _1 V 


n.a. 


n 


reduced shear: 11 = x~ 2 H 


n.a. 



TABLE I: Symbols and their meanings. 

dynamics is governed by the Klein-Gordon equation. For 
the background, it reads: 

(p = -2H0-a 2 V'{ip), (A31) 

while the perturbation obeys the equation of motion 

X = - 3$) - 2a 2 V'((p)^ 

- {a 2 V"{Lp) + k 2 )X - 2HX. (A32) 

From the energy momentum tensor for the 
quintessence field 

T» v = <p«<p iV - 5» (j<p' a <p, a + V(<pfj > (A33) 

using ip — (p + X and the longitudinal gauge metric, one 
gets 



ST 



(ion.) 




£ 2 $ _ Xtp ) - V'(ip)X Q (A34) 



5T\ {lon - ] = -a- 2 k(pXQ l 



(A35) 



Using the definition of A, equation l|A22|) in longitudinal 
gauge and p q +p q = cT 2 (p 2 one can read off from equation 
the gauge-invariant expression 



evolution equations 



(1 



2a 2 V'(<p) 
6aaV'((p) 



A„ 



1 + Wn 



3$ 



kip 



k V„ 



(A38) 



w a 



and 



V q = k 



3$ + * 



Wn 



2HV a 



(A39) 



Equation l|A38[) depends on the specific quintessence 
model through V and (p. We can however make progress 
in the case of nearly constant w q : Many quintessence 
models have solutions for which ip approaches an attrac- 
tor solution irrespectively of its initial value. For these 
tracking quintessence models 0, 0, 0] , the equation of 
state of the quintessence field w q is nearly constant dur- 
ing radiation domination. We will use this vanishing of 
w q in the following to derive relations to simplify equa- 



tion (IA38 



Considering a 2 ip 2 



(1 



w q )p v 



using the Friedman equation 3a MpTC 
cp= [3(l+w,)fi,]5MpH, 

and hence 

(£ _ _d_ A _ If!, H 

~ dr ULf ~ 20, + H' 

where we have neglected a term involving w q . We will 
in the following assume that at early times, the universe 
expands as if radiation dominated. In this case, Ti. = t _1 
and inserting the above equation (|A41|) into the equation 
of motion (j A3 1(1 , one finds 



it follows 
p that 

(A40) 
(A41) 



a 2 V _ 3(1 -w„) 
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2t 



(A42) 



Using this relation (|A42|) . the evolution equation for A q 
becomes 



3(w q 



A 
3- 



., -3(1 + 0$ 
2 



x 



(1 + W q )V q 



(A43) 



3K-1) 

whereas the one for the velocity remains almost unaltered 
while we move to the reduced velocity V q : 



V a 



A„ 



l + W 



3$ + * 



r- x V a 



(A44) 



A q = (1+Wg) 



3$ - * + X0~ 



+ XV'((p)p q 1 '. (A36) Note that T q does not usually vanish. Instead, we obtain 



In the same manner, one gets from equation l|A35|l and 
the fact that v ( io "f ) = V the relation 



w q T q = (l-c 2 s{q) ) 



A q - 3(1 + w q )<P + 3-(l + w q )^- 



V q = k(p X. 



Taking the time derivative of equations i|A36(l and (|A37|) 
and using the equation of motion i|A32Jl , one obtains the 



(A37) with the sound speed of quintessence given by 



la w„ 



Wn 



(A45) 



(A46) 
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Matter and Radiation 



Setting w = c 2 s = T = in equations (|A29|I and HA30|> . 
we obtain the cold dark matter evolution equations 



U 



-kxV r . 



(A47) 
(A48) 



The multipole expansion of the neutrino distribution 
function can be truncated beyond the quadrupole 
at early times. In terms of density, velocity and shear, it 
is given by [3(], H|J 



A, 



--kxV u , 
k f 1 A 



6 



21L 



(A49) 
* <i> j , (A50) 
(A51) 



Deep in the radiation dominated era, for which the initial 
conditions here are derived, Compton scattering tightly 
couples photons and baryons |2y,|33- The coupling leads 
to Vh = Vy and the evolution equations become [2(| 



A 7 = 


4 - 
— — kxVy, 


(A52) 








U 7 = 




(A53) 


A b = 


— k X Vry. 


(A54) 



As the photon quadrupole and all higher photon multi- 
poles are suppressed during tight coupling, it follows that 
$ is given from Einstein's equation by 



$ = -tf-fi„n„, 



(A55) 



where we have used the Friedmann equation. Finally, the 
Poisson equation i|A26(l in terms of the various species is 



E n a (Aa + 3(i + w a )v c 

a=c,f>,7,f ,q 



E 3(l + w a )Sl a + *f 



21 <!,ll,. (A56) 



where the index a runs over all species. Rewriting the 
evolution equations i|A47(l - (|A54|) in terms of d/dlna; 
and replacing <I> by means of dsni, one arrives at J2J- 

CD- 



APPENDIX B: EARLY TIME QUINTESSENCE 
FIELD FLUCTUATIONS 



While throughout this work, we describe quintessence 
perturbations by the variables {A qi V q }, one could in- 
stead use the field fluctuation and its time derivative 



Quantity 



Scaling behaviour 



9 
V 
V" 

\ adiab. 
-^adiab. 



oc r 



-(l+3ui„)/2 
-(7+3tu 9 )/2 
-4 



oc r 
const. 

x T (l-3w q )/2 



TABLE II: Tracking quintessence in the radiation era: Scaling 
handbook. 



{X,X}. In this section, we will give analytic expres- 
sions for X and X in the case of tracking quintessence 
for super-horizon modes. We will do so assuming that \& 
and <& are at least almost constant. As this is not the case 
for CDM isocurvature and baryon isocurvature, the fol- 
lowing steps do not apply in these modes. Furthermore, 
we will assume that the universe expands as if radia- 
tion dominated during the time of interest. In this case, 
H = r -1 , n q oc r't-^wq anc j hence by means of equa- 
tion 1A40|I <p oc T^^ < - 1+3Wq \ Using this, we infer from 
equation (|A42|) that V cx T ^^ 7+3w i'> . In addition, a 
straightforward calculation using (|A41|) and (|A42(I yields 



, rr _dy^dr = h 1 _ Wq) ( 7 + 3Wq) _ (B1) 



dr dtp 4 



The equation of motion for X (|A32() contains a term 
(p (jfr — 3$^ , which by assumption we may drop. In ad- 
dition, we see from equation I|IJ1|I . that for super- horizon 
modes, a 2 V" ^> k 2 (except for w q very close to 1), and 
hence the equation of motion reduces to 



X = -2a 2 V'® - a 2 V"X ~ 2-X. 

a 



(B2) 



Using the power law behaviour in r of V, V" and a, 
as well as equations (IA42II IjBlll . one hnds the particular 
solution 



(B3) 



as well as two complementary solutions that may be 
added to obtain the general solution 

X(T) = I* <p + C X T -*(Wl-4aW) 

+ C2T -|(l+Vl-4^r^"). (B4) 

The mode proportional to ci is at least as rapidly decay- 
ing as the one proportional to c\ . Using the explicit form 
of Aa 2 T 2 V", equation JETJ|, we find that y/l - Aa 2 T 2 V" 
is imaginary if w q € [— 1(1 + y/E),— 1(1 — v6)], which 
holds for all scalar quintessence models of current inter- 
est. Hence, the complementary modes decay oc \j\fr in 
an oscillating manner. 
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Coming back to the dominating particular solution 
(|lj3[l , Figure shows that the accuracy of this analytic 
result is indeed high at early times, when compared to 
numerical simulations. 

Inserting the solution i|B3l) and its time derivative into 
equation IJA36J) . we find the simple expression 



While the conclusion regarding Uq and U\ are still the 
same as in section l|ll B[) . we see that quintessence may 
introduce a correction 



A q = 3(1 + w q ) 



2 



(B5) 



which is just a restatement of eqn. 123( 1 and (|24() . Hence, 
the energy density contrast in tracking quintessence mod- 
els remains constant on super horizon scales, provided the 
gravitational potentials are constant to good approxima- 
tion. 



U a 



[A 



iwq)!]- 1 A q U . 



(C6) 



APPENDIX C: EXTENDED MATRIX 
FORMULATION 

For simplicity, we have limited the discussion in section 
(|IIBI| to cases where either w q = 1/3, or quintessence 
contributions to A(x) are neglected. Here, we will discuss 
cases for which w q < 1/3, while the background expands 
radiation dominated. In this case, flq cx rf 1 " 3 ^) and 
we can split the matrix in three parts according to their 
scaling with x: 

A{x) = A + x + x^- 3w ^ A q . (CI) 

Again, Equation (Q will lead to a solution vector of the 
form 

U(x) =U a + xUi+x {1 - Zw " ) U q . (C2) 

Substituting this into Equation Q and keeping only 
leading orders in x, we get 

A Uo = 0, (C3) 
AxUo + AoUx = U u (C4) 
AqUo + AoUq = (l-3io,) 17,. (C5) 



This contribution x^ 1 ^^^ U q could in principle domi- 
nate over xU\ for w q > 0, fl q > fl c . However, the con- 
tribution is only of interest for the CDM isocurvature and 
baryon isocurvature modes, as it is otherwise negligible 
compared to the constant order. Yet for CDM isocur- 
vature and baryon isocurvature, A q Uq = 0. Therefore, 
the discussion below applies, leading to U q = for CDM 
isocurvature and baryon isocurvature modes. One order 
higher in x, there may be a contribution. Yet this is 
in any case a higher order contribution, which we may 
neglect. 

Finally, we briefly discuss the case of vanishing Uq. 
This only concerns possible subdominant modes. Equa- 
tion (|C4|) then yields Aq U\ = U\, i.e. U\ is an eigen- 
vector of Aq with eigenvalue A = 1. As Ao does not 
have such an eigenvector, we are led to conclude that 
Equation does not have a regular solution involving 
Ui, if Uq = 0. Turning to Equation (|C5|) . we similarly 
conclude that U q needs to be a eigenvector of Aq with 
A = (1 — 3w q ) for vanishing Uq. For w q < 1/3 this is 
once again excluded and for w q — 1/3, we just regain the 
results of section Hi Bl 
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